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Abstract. Let H be the discrete Schrodinger operator 
Hu(n) := u(n — 1) + u(n + 1) + v{n)u(n), u{0) = acting on 
P(Z"'") where the potential v is real- valued and v{n) ^ as n — > oo. 
Let P be the orthogonal projection onto a closed linear subspace 
C C P(Z+). In a recent paper E.B. Davies defines the second order 
spectrum Spec2(iJ, £) of H relative to C as the set of z S C such 
that the restriction to £ of the operator P{H — z^P is not invert- 
ible within the space £. The purpose of this article is to investigate 
properties of Spec2(-ff, £) when C is large but finite dimensional. 
We explore in particular the connection between this set and the 
spectrum of H. Our main result provides sharp bounds in terms 
of the potential v for the asymptotic behaviour of Spec2(i?, >C) as 
C increases towards 1^(Z+). 

1. Introduction 

Let the discrete Schrodinger operator H be defined by 

Hu{n) := u{n + 1) + u{n — 1) + f (n)-u(n), m(0) = 

acting on 1^(Z+) where v : Z+ — > R. In j2], E.B. Davies investi- 
gated the concept of resonance of H associated to a closed subspace 
C C Let P be the orthogonal projection onto C He called the 

isolated points of the second order spectrum of H 

Spec2 {H; £) := {A G C : P{H - \fP\C is not invertible}, (1) 

resonances of H relative to C This definition is motivated by graphical 
and numerical connections between other notions of resonance adopted 
to this context, and isolated points of Spec2 [H] C) (cf. [3 section 9]). 
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This was observed by Davies in the case when the dimension of L is 
low (< 30) and v has finite support. 

In Davies also pointed out that Spec2 C) contains information 
about the spectrum Spec (if) of H. In jH] E. Shargorodsky developed 
this idea for general linear operators T. This began the study of sec- 
ond order spectra as a projection method for localizing the spectrum. 
Among other interesting results on geometrical properties of Speca, 
Shargorodsky showed that for any bounded self-adjoint operator T, 

I J lim Spec2 (T; £fc) n M = Spec (T) (2) 

where the union is taken over the set of all sequences (Ck) of subspaces 
of the domain of T such that the orthogonal projection onto Ck 
converges strongly to the identity, see theorem 21]. Hence the second 
order spectra of H might be useful for approximate computation of 
Spec (H) when v is bounded. 

The purpose of the present paper is to continue investigating the 
connection between the spectrum and second order spectra of H. We 
present a detailed description of 

Spec2 {H; k) := Spec2 {H; Ck) 

for large k when 

A := span{5i, . . . ,4}, : = 

and there exist constants a, r > such that 

\v{n)\<^ n = l,2,... (3) 

As V is real-valued and bounded, if is a bounded self-adjoint operator 
so we are in the situation covered by Shargorodsky's results. Since © 
holds for finite rank v, our discussion includes the potentials investi- 
gated by Davies in [2]. 

By construction, Ck is fc-dimensional. Then the quadratic operator 
pencil in is a X matrix and 

Spec2 {H; k) ={XeC : det[Pk{H - XfPklCk] = 0}. 

It is easy to see that the determinant at the right hand side is a polyno- 
mial in A of degree at most 2k so that Spec2 {H; k) consists exclusively 
of at most 2k relative resonances of H. 

Since our interest is in how the Spec2 (ii; k) are related to the spec- 
trum of H, we describe the latter. The classical theory allows us to 
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find Spec(iJ) under the hypothesis © quite efficiently. For this we 
decompose 

H = Ho + V 

where i^o is H without potential and Vu{n) := v{n)u{n). The spec- 
trum of Hq is pure absolutely continuous and equal to the interval 
[-2,2]. Since 

oo oo 

Tr(V) = Y^\{VS^,6^)\=Y,\v{n)\ 

n=l n=l 

oo 

n=l 

\^ is a trace class operator. Then, by virtue of Pearson's theorem, 

Spec,,, (H) = Spec,, (H) = [-2, 2] 

so Spec (if) consists of the interval [—2, 2] together with a discrete set 
of isolated eigenvalues of finite multiplicity which can only accumulate 
at ±2. A recent adaptation by F. Luef and G. Teschl of the Sturm- 
Liouville theorem to the discrete context, cf. shows that in fact if 
r > 1 in 0, then H has only a finite number of eigenvalues. 

The results we will discuss in the forthcoming sections are motivated 
by the following proposition. 

Proposition 1. Let V = V* be a compact operator acting on 1^(Z+). 
Then for all z ^ Spec {Hq + V) there exists k > such that 

z ^ Spec2 (Hq + V;k) , k>k. 

In other words, roughly speaking 

S := lim Spec2 {H; k) C Spec(i7) . (4) 

A;— >oo 

Notice that we are less restrictive on the requirements for V . In com- 
paring with this result says that our particular choice of Ck does 
not lead to spurious points in approximating Spec (if) . In section|21we 
will show proposition Q using an argument that will be crucial for our 
latter work. It consists in considering Speca (if; k) as the perturbed 
spectrum of a convenient Toeplitz matrix. 

Our main results are theorem|3]and corollary |S1 of section^] There we 
give a precise meaning to the limit in (0} and provide sharp estimates 
on the rate of convergence in terms of v when it satisfies 0. 

The proof of theorem |3] requires to develop some results on the sta- 
bility of Toeplitz operators with smooth symbols. This is carried out in 
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section 01 It will become clear from the beginning that these techniques 
can be extended to more general situations, hence we give a general 
treatment to the topic. Section El is closely related to chapter 2 of the 
monograph 1 on large Toeplitz matrices. Since we failed to find a 
reference for theorem |21 and corollary El we include proofs of both. 

We complete our discussion by considering in detail some numerical 
examples in sectional These are closely related to those of [3 section 9]. 



2. The finite section method 

The finite section method for perturbed Toeplitz operators was de- 
veloped by I. Gohberg, LA. Feldman and H. Widom in the 1970's (see 
|3] and jU]). Through this paper we follow closely their ideas, in par- 
ticular we follow the notation and procedures of the introductory text 
on the subject by A. Botcher and B. Silbermann, PP. This section 
is devoted to discuss the basic notation and results in the Gohberg- 
Feldman- Widom approach which will be the core of our latter work. 
This will lead us to a proof of proposition HI 

Below and elsewhere we identify the truncation 

PkTPk\Ck '■ Ck — ^ C,k 

of any linear operator T : 1^(Z+) — > 1^(Z+), with its matrix repre- 
sentation acting on C*=. We will often denote := PkTPk\Ck. The 
following definition is taken from jH chapter 2]: we say that the se- 
quence {PfcTPfc} is stable (with constants ko and M) if and only if, 
there exists ko > such that Tk is invertible for all k > ko and 

sup ||T;ri|| < M < oo. 

k>ko 

Notice that if {Pk{H — z)'^Pk} is stable then z E C falls outside 
Spec2 {H; k) for large k. 

If T is a Toeplitz operator with continuous symbol, then {P^TPk] is 
stable if and only if T is invertible (cf. ^1 theorem 2.11]). Furthermore 
(see corollary El below or [H theorem 2.16]), if {PkTPk} is stable, K is 
compact and T + i^' is invertible, then {PkiT + K)Pk} is stable. 

The following argument is crucial to out later work. Decompose 



{H - zf = {Ho - zf + {Ho - z)V + V{Ho -z) + V^ 
= T{(5) + K{z) 
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where T{(3) is the Toephtz operator whose symbol is 
f3{t; z) := - 2zt-^ + (2 + z'^) - 2zt + 

= {r^ - zr^ + -zt + i) 

= {t-'-z + tf, 
teT:= {e*'^ : -vr < ^9 < vr} 

and 

K{z) := {Ho - z)V + V{Ho -z) + V^- \5r){5^\ 
= HqV + VHo + V^- 2zV - \6i){6i\ 

is compact. According to the above, the stabihty of {PkT{P)Pk} im- 
phes stabihty of {Pk{H — z)'^Pk} when z ^ Spec (if) . This is the key 
argument in the proof of proposition ^ 
Proof [of proposition The symbol 

(3{e'^;z) = {2cos^ - zY, -tt < ?9 < tt. 

Then /3(T; z) is an open curve without loops so therefore, by virtue of 
Gohberg's theorem on the spectrum of Toeplitz operators with contin- 
uous symbol (see for instance [H theorem 1.17]), SpecT(/3) = P(T;z). 
Thus T{j3) is invertible if and only if 2: ^ [—2,2]. For all such z, 
{PfcT(/5)Pfc} is stable. Clearly [H — z^ is invertible if and only if 
z ^ Spec (if) . Hence, 

{Pk[T{f3) + K{z)]Pk} = {Pk{H - zfPk] 

is stable for all z ^ Spec (if) . ■ 

3. Stability of Toeplitz operators with smooth symbol 

In this section T{a) denotes the Toeplitz operator associated to the 
symbol a : T — > C. If g > 0, stands for the set of continuous a 
such that 

|n|'^|Q'(n)| < 00 

where a{n) denote here and elsewhere the (non-normalized) Fourier 
coefficients of a. It is well known that is a Banach algebra with 
pointwise algebraic operations and norm ||a|| := + 
According to Wiener's theorem (cf. [S]), if a G and ^ tt(T), then 

Throughout this section we will assume that 
i) The symbol a G for some q> 1. 
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ii) The sequence {PkT{a)Pk} is stable with constants denoted by 
ko and M: that is Tk{a) is invertible for all k > ko and 

sup,>,„ \\T^\a)\\<M. 

The continuity of a and stability of {PkT{a)Pk} guarantee that T{a) 
is invertible. Under these conditions it is well known (cf. pp.32- 
34]) that TfT^ia) converges in the strong operator topology to T~^{a). 
Therefore it is legitimate to estimate the columns of the matrix of 
T~^{a) from the columns of T^^(a) for large k. The following result 
states that if we not only assume continuity of the symbol but the 
stronger condition i), the error in the estimation of the first k/2 columns 
of T-i(a) from T^^{a) is 0{k'i-'). 

Theorem 2. Let T{a) satisfy Then for all 1 < p < q there 

exists a constant c > independent of k, such that 

\\T^\a)5,-T-\a)5,\\<^^ 

provided k > ko and 1 < j < k/2. 

Below we will compute an explicit upper bound for c in terms of a, 
M and \\T-\a)\\ =: M. 

This theorem improves theorem 2.15]. It allows us to show a ver- 
sion of the result mention in section |21on stability of compact perturba- 
tion of stable sequences, when the perturbed operator T is Toeplitz, its 
symbol satisfies i) and the perturbation is a trace class band matrix. 
This is the content of corollary El Let 

K,,i = {K6i,6,) I, J = 1,2,... 

Kj^i = / < or j < 0. 

We assume that K satisfies the following hypotheses, 

a) There exists L > such that Kj^i = for |j — /| > L. 

b) There exist 6 > and r > 0, such that 



b 

f 



Kj,i<-^, I, J = 1,2, 



Corollary 3. Let T{a) satisfy i)-ii) and let K satisfy a)-b). Fix 
I < p < q, let p := min{r, (2p — l)/2} and let c be as in theorem\^ If 
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T{a) + K is invertible, then Pk{T{a) + K)Pk\Ck is invertible provided 
k > 



4:L 



max 



cV3\\K\\ + hT-^{M + M) \\{T{a) + K)-^T{a) 



(6) 



Proof . Throughout the proof we assume 1 <p < q and k > kg. By 
virtue of ii) , 

I + T-\a)K = T~\a){T{a) + K) 

is invertible. Let 

l/e:= \\{I + T-\a)K)-^\\ = \ma) + K)-^T{a)\\. 
Then for all m e 12(Z+), 

e\\Pku\\ < \\{I + T-\a)K)Pku\\ 

<\\{I -rT-\a)P,K)P,u\\ + \\T^\a)P,K -T-\a)K\\ \\P,u\\. 

The operator K acts on as follows, 

oo L 

Ku^^ ^ Kn,n+mU{n + m)(5„. 
n=l m=—L 

Let Rk := T^\a)Pk - T-\a). Then 

oo L 



\\RkKu\\ 



n=l m=—L 

L cc' 



,n+m I \\RkSn\ 

m=—L n=\ 

For each m = —L, ...,!/, 

oo 

,n+m 



n=l 



1/2 



1/2 



< [J2Hn + m)\^ [J2\Kn,n+m\'\\RA 



,n=l 



vn=l 



k/2 



1/2 



^ Ikll ( y^jKn,n+m\^\\RkK\\'^ + ^ |-f^n,n+mH|-Rfc(^n| 

n=fe/2+l 



n=l 
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By virtue of theorem |2l 



fe/2 „2sr^k/2 

E\K A l|2 ^ l^n=l K^", 



n+m I 
n=l 

12 -1 7|| 7^112 



By virtue of a)-b) and by the definition of Rk, 



n=A;/2+l n=fe/2+l 



n+m I 



n=k/2+l 



< 



2 



^ 6^2^-(M + M)^Er=./2+il/^ 



< 



6222'-(M + M)27r 



6A; 



2r 



Then, if k satisfies (P), 

For all such k, 

{6/2)\\Pkv\\ < \\{I + T,\a)PkK)P,v\\ 

<\\T,\a)\\ WiUa) + PkK)P,v\\ 
= \\T,\a)\\ \\Pk{Tia) + K)Pkv\\ 
<M\\PkiT{a) + K)Pkv\\, 

so that Pk{T{a) + K)Pk\Ck is injective and thus invertible. Notice that 
dim Ck is finite. ■ 
In section 0] we will employ this corollary to estimate sharp bounds 
for the error in the limit (j3)). 

We now prove theorem |21 For this purpose we use the following 
convention: Qk := I — Pk, 



Wku(n) :-- 



u{k — n + 1) 1 < n < k 
n > k 
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a := a{ir^\ H{a) is the Hankel operator generated by a (cf. P-13]), 

A{a) := T-\a) - T{a-^) 

and 

B, := P,T-\a)Pk + WkAia)Wk. 

Assumption ii) and the fact that a is continuous, ensure that A (a) is 
well defined. The following results are due to Widom, but a proof can 
be found in pp. 40-42]: If a E L°°{T) and T{a) is invertible, then 

1) A(a) = T-\a)H{a)H{a-^) = H{a-^)H{&)T'\&), 

2) Tk{a)Bk = Pk- PkT{a)QkA{a)P, - WkT{a)QkA{a)Wk. 

Proof [of theorem Throughout the proof we assume that 
k > kQ and 1 < j < k/2. Decompose 

T,-\a)=Bk + Ck, 
Ck := T^\a) -Bk = T-\a){Pk - Tk{a)Bk). 

Then 

T,-i(a)5, - T-\a)5, = [T,-\a) - P,T-\a)P,]6, + [PkT-\a)Pk - T-\a)]S, 

= [Bk + Cfc - PkT-\a)Pk]6, + [PkT-\a)Pk - T-\a)]6, 

= Ck6j + WkA{a)WkSj + [PkT~\a)Pk - T-\a)]5j. 

(7) 

Fix 1 < p < q. In order to find the parameter c, we will estimate the 
norm of the three terms in the sum at the end of ((Tj). For this, let 

Ci := |n|^|Q;~^(n)|, 
C2 := [XII^I'I^^HI') ■ 

VneZ / 

Since a G and ^ tt(T), both Ci and hence C2 are finite. By 
definition of Qk and H{a~^), 

\\Q,H{a-^)\\<\\a-\.)- 5^He*""l|oo 

n<k—l 

\n\>k 

^ EH>fcl^l^l«~H^)l ^ (h 

~ kP ~ kP' 
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Hence by virtue of 1), 



lin A^^^ll <r ^i\\H{&)T-\a)\\ c,M\\a\ 

\\QkHoi)\\ < < 



kP ~ kP 

In a similar manner one can show 

CiM||a||oo CiM||q;||oc 



\\QkHa)\\ < 



kP kP 
We estimate the norm of Ck- By virtue of 2), 

Ck^-T^\a){n{a)Bk-Pk) 

= T,\a){PkT{a)QkA{a)Pk + WkT{a)QkA{&)Wk). 



Thus 



llCfcll < \\T,\a)\\ \\PkT{a)QkA{a)Pk + WkT{&)QkA{a)Wk\\ 
< M (||a|U||gfcA(«)|| + ||a|U||gfeA(a)||) 
2MciM||a;|'2 



loo 



< 

kP 

The second term is 

\\WkA{a)WkSj\\ < ||A(a)4-,+i|| 

= \\T-\&)H{d)H{a-')6k-j+i\ 
< M\\a\U\H{a-')Sk^j+i\\ 

. 1/2 

\n>k—j 

1/2 

<M||a||oo ( Yl I^H''^ 

\n>fe/2 



\a-HnV^ 



< 



kP 

2PM\\a\\ooC2 



1/2 



< 

~ kP 
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The third term is 

\\[PkT-\a)Pk-T-\a)]5,\\ = \\PkT-\a)5, - T-\a)6,\\ 

= \\QuT'\a)6,\\ 

= \\Qu[Ka)+na-')]5,\\ 

< \\Qu^{a)5,\\ + \\QkT{a-^)5,\\ 

Ci\\a\\ooM 
'kP 



< n; ^ \\QkT{a )dj\ 



1/2 



kP 



ci a ooM I \- — 
^ TTp + 2^ I" 'Wl 



n>k—j-\-l 

1/2 

2 



.n>k/2 



< 



Ci||a|LM + 2Pc2 



n>k/2 



kP 

(11) 

Hence the conclusion of theorem |21 can be recovered from (0), (0), (fTUj) 
and (fTT|) by putting 

c = M||«||oo(2Mci||a||oo + 2^02 + ci) + 2^02- ■ 



4. The second order spectra of H 

We are now ready to state and prove the main resuhs of this paper. 
These roughly say that for all k large enough, Spec2 {H; k) is contained 
in a small neighbourhood of Spec (if) with diameter of order some 
negative powers of k. As we previously mentioned, applying corollaryOl 
to the decomposition (0) will play a crucial part into the proofs. 

Theorem 4. // (jSj) holds, then there exists a constant b{r) > inde- 
pendent of z, such that 

Spec2{H;k) c{zeC : dist (z. Spec if) < 6(r)/A;2r/(2r+35)|^ 
for all k. 



This automatically implies: 
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Corollary 5. Suppose that r > can be chosen arbitrarily large for 
f m 0. Then for all < q < 1 there exists a constant b{q) > 
independent of z, such that 

Spec2 {H;k) c{zeC : dist (z, Spec i/) <b{q)/k''}, 

for all k. 

The rest of this section is devoted to proving theorem 01 We first 
introduce some further notation. Let 

d{z) := dist {z, [—2, 2]) and d{z) := dist {z, Spec H) . 

Then d{z) > d{z). For z ^ [—2,2], we denote by T{P~^) the Toeplitz 
operator whose symbol is 

p-\t; z) = {t-^ - 2zt-^ + (2 + z^) - 2zt + t^Y^. 

Since [3^^{T;z) is an open curve and ^ /?~^(T;2;), then T(/5~^) is 
invertible. Put Ms := \\T-^{(3-^)\\. 

We break the proof into various steps. The aim of these steps is to 
show that for all i? > 0, there exists 6 > independent of z, such that 
Pk{H — z^Pk is invertible provided 

Spec (if), d{z)<R 

and 

k > ^ 

- [J(2)](2r+35)/2r- 

Since the Specg [H; k) are bounded uniformly in k (cf. [3 theorem 18]), 
this assertion implies theorem 0] Below we assume without further 
mention that z Spec (if) and d{z) < R. The various constants bj 
that appear below are independent of z but they may depend on R, p 
or V. 

Step 1: we estimate M and M2 for the symbol under consideration. 
When z ^ [-2,2], 



s± 2 

are such that C+ = C-^ and \C±\ 7^ 1. One of these numbers is inside 
the unit disk whereas the other is outside it. Denote by ( the one that 
is outside and put 

P±{t;z) :=l(t±i-c)^ 

Then 

(3{t;z)=(3+{t; z)(3.{t;z), 
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can be extended analytically inside the unit disk and can be 
extended analytically outside it. This is the so called Wiener-Hopf 
factorization of (3, so that (see for instance jll theorem 1.15]) 

Hence 



M = \\T-^' 

<IIW;^)II Wnpz') 



^+ 1 1 oo 1 1 P— 1 1 oo 



< 



icrsup|t-cr'sup|i/t-ci 
icrsup|t-cr' 

bo 

d{zY' 



-2 



By writing the corresponding Wiener-Hopf factorization for /5 ^ we 
obtain in a similar manner 



M2<||T(/3+)|| ||T(/3_: 



sup|t-Crsup|l/t-CP < hi. 



Is I teT teT 

Step 2: we now want to estimate the constants Ci and C2 in the proof 
of theorem 121 in terms of z. Notice that since /? is a trigonometric 
polynomial, P and hence P~^ belong to for all p = 1, 2, . . .. It is 
easy to show that 



dP 
d3P 



p+2 



, (2cos^9-;z)^' 

r=3 ^ ' 

where the (pri'd) are smooth, bounded and independent of z. Then 

1/2 / X 1/2 



C2 



< 



< ho 




< 



p+2 _^ 



< 



(2 cost? - z)'' 

&4 



1/2 

d?9 



1/2 



r=3 
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Similarly 

1/2 / \ 1/2 

/ ^ 1 \ 

< 



< ^ 



d{z] 



p+3' 



Step 3: estimation of M. This only makes sense for k large enough. 

Lemma 6. For all R > and 1 < q < 2, there exists constants 
Bq, bj > independent of z, such that Tk{f3) is invertible and 

\K\m<^,, (12) 

provided 

z^[-2,2], d{z)<R 

and 

b7 

Proof . By virtue of lemma 2.9], we know that Tk{(3) is invertible 
if and only if 

QfcT-i(/5)Qfc|RanQfc 
is invertible, and in this case 

Ti:\m = 

PkT-\(3)Pk - PkT~\p)Qk{QkT~\P)Qk\R^nQk)-'QkT-\P)Pk. 
The truncation 

(13) 

where 

A{(3)=T-\f])-T{(3-') 

as in section|21 The matrix of QkT{i3~^)Qk\RsinQk is the same matrix 
T{j3~^). Then, since ^ /5~^(T), the former is invertible and 

{QkT{(3-')Qk\RanQky' 

has the same matrix as T~^{f3~^). Thus 

\\{QkT{(3-')Qk\Ra.nQk)-'\\ = \\T-\(3-')\\. 
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Let p = 1, 2, . . .. By virtue of 

\\Q.mQ.\\<"-^^^^ k>i. 

According to the steps 1 and 2, 
Then, if 

k > ^ 



\\{QkT{(3-')Qk\RanQi,y'QkA{l3)Qk\\ 

< \\{QkT{r')Qk\R^nQk)-'\\ ||QfcA(/3)Qfc| 

\T-\p-')\\c,M\\(3U 



< 



kP 

CiM||/3|UM2 



< 1/2. 



kP 

By virtue of (fT^. for all such k, Q ^T'^ {P)Q k\Ra.n Q ^ and so Tfc(/3) are 
invertible, and 

\\T,-\m<\\T-\m+nT-\f3)rM2 

r2f. ^ bo , 2bobi 



<M + 2M^M2 < -TTJ + 



d{zY d{zf 

Take p large enough to complete the proof. ■ 
Below M denotes the right hand side (fT^. 

Step 4: We are now ready to complete the proof of theorem For 
this we use corollary El According to (jSj), 

Pu{H - zfPk = + PuK{z)Pk, 

where 

K{z) = HoV + VHo + V^ -2zV -\6i){6i\. 
It is easy to see that the matrix entries of K{z) are 

{K{z)6n, 6n) = 6i{n) + v\n) - 2zv{n) 
{K{z)6n, Sn+i) = v{n) + v{n + 1) 
{K{z)5n+i, 5n) = v{n) + v{n + 1) 

{K{z)6i,6m) = if|/-m|>2. 

By virtue of 0, this ensures that K{z) satisfies conditions a)-b) of 
section El 
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Fix p = (2r + l)/2. According to steps 1-3, the constant c of theo- 
rem |21 for this particular symbol is 

c = M||/?|U(2Mci||/?|U + 2^C2 + ci) + 2^C2 



d{zy \d{z)P+^^ d{z)P+^ d{z)P+^ J d{z) 



p+2 



< ^15 ^ _ ^15 



Also \\K{z)\\ < 6i6 and 

mp)+K{z))-'Tm\<\\{H-z)-Y 

^ i2+\z\r ^ by, 



d{zf ~ d{zY 
Then 

[cV?>\\K\\+b2'^{M + M)]\\{T{p) + KY^T{p)\\ < 

bi8 bi8 



^(^)(2p+34)/2 ^(^)(2r+35)/2 

Thus, by virtue of corollary El and lemma IHl there exists 6 > such 
that Pk{H — z^Pk is invertible provided 

z0 Spec (if), d{z)<R 

and 

b 

k > 



[^(^)](2r+35)/2r 

This completes the proof of theorem |3] ■ 

In principle one can follow track of the constants bo, . . . ,bis for par- 
ticular potentials. Nonetheless one should be aware that b increases its 
value faster than 2*" as r increases towards infinity. 

5. Numerical examples 

In this final section we discuss corollaryElfor rank one potentials from 
the numerical point of view. In particular we show some computations 
of second order spectrum for this kind of potentials and compare these 
with Spec (if) which in this case can be found explicitly. The results 
below extend those of \^ example 20]. 

Everywhere in this section, we assume that v is such that 

'^^^^ I otherwise, ^^^^ 
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for some j = 3, 4, . . . and a > 2. In the proof of the following result we 
use some properties of the transfer matrix associated to the difference 
equation 

Hu{n) = \u{n), u{0) = 0, u{l) = 1. (15) 

See |3] and [7j for some recent results on the relationship between spec- 
tral properties of H and the transfer matrix for slow decaying poten- 
tials. 

Proposition 7. Let v be as in (jl4j) . Then the discrete spectrum of H 
consists of an isolated eigenvalue Xa, such that 

a < \a < a + 2/a. 

Proof . Since a is positive, any eigenvalue Xa of H should be Xa > 2. 
Let A > 2 be a solution of the eigenvalues equation (jl5|) where u{n) is 
a sequence with (at the moment) no constraint of growth at infinity. 
For all n > 1, 

u{n) \ _ f 1 ^ fu{n - i; 

u(n + l)l I— 1 X — v(n)l \ uin) 



T 



1 
-1 A 



Let 



Then for all n ^ j, 



/u{n +1)/ V u(n) 
The Jordan decomposition for T yields 

T = u(^^ ]u-\ 



u(n) \ ^fu{n-l), ^^g^ 



where fi±{X) = fi± := (A ± VA^ -4)/2, 

U=(^ M and U-^ = (fi- - ( 

By applying inductively (fT6|l for n from 1 up to the step j — I, we 
obtain 

u{j) ) - ^ \u{i)) - ^ ^ ^rV VI 
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Then 

' u{j) A / I \fu{j- 1) 
Mj + l)J l-l X-aJ\ u{j) 



3 3 



Hence, for all n > j 

+ Wi + 1) 



n— 7 



where 



- + 1) , 
wilX) = wi :— and 

W2(A) = W2 := ^ —. 

If M is an eigenvector in then necessarily u{n) — > as n — >• oo. 

Hence 

wi — > and ii"L~''w2 
as n — > oo. Since A > 2, then /i+ > 1 and < /i_ < 1. Notice that wi 
and ■u;2 are independent of n. Thus in order for u e P(Z+), necessarily 
wi = 0. 

By substituting u{j) and M(j + 1) in the identity for wi, 

^ l^-il^- - /^+) - (A - a)(/xi - + (/^r^ - 1^^+^) 

Notice also that = /x^. Then, a straightforward computation 

shows that wi — 0, if and only if q{n+) — where 

q{x) := x^^^^ - ax^^ - x^^'^ + a. 

It is easy to see that q'{x) vanishes, if and only if x = or 



_ 2ja ± V(2ja)2 + 4(2j + l)(2j^ 
Aj + 2 

Furthermore 

1) For all j > 3, a;_ < and x+ > 1, 

2) q{0) = a > 0, 
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a 


j 


K 


Estimation of using Speca {H] 60) 


3 


3 


3.60362098809610 


3.60362098809609 - 0.00000008894009? 


3 


6 


3.60554979388607 


3.60554979388608 - 0.00000005798745z 


3 


9 


3.60555127432257 


3.60555134832088 


3 


12 


3.60555127546311 


3.60555133471164 


3 


5 


3.60553511316735 


3.60553516591671 


6 


5 


6.32455524824927 


6.32455533827497 


9 


5 


9.21954445506085 


9.21954445506086 - 0.00000009147845z 


12 


5 


12.16552506041796 


12.16552524742606 



Table 1. Aa for selected values of a and k. 



3) q{l) = and 

4) q'{x) < for < X < x+. 

Hence, it is not difficult to see that q{x) has only one root Xa in the 
interval (1, oo). Since q{a) = a — a^^~^ < and 

q{a + l/a) = > 0, 

a + a-^ 

necessarily a < Xa < a + 1/a. Now, 

/i+(A) = Xa, 

if and only if 

\ = \^-= Xa+ l/Xa- 

Therefore Aa is the only possible eigenvalue of H and 

a < Xa < a + 2/a. 

Finally we must show that A^ is indeed an eigenvalue of H. By virtue 
of dm), for all n> j 

U{n) = H":~\\a)w2{\a)- 

Since W2{Xa) 7^ is fixed in n and < /i_ < 1, m is an P(Z''") eigenvector 
of H as we required. ■ 

In the third column of table Q we show the value of A^ for selected 
a and j. We found the roots of the polynomial q{x) using the internal 
algorithm "roots" that the package Matlab provides. More complete 
data can be found without difficulty by this method. 

In figures El and 01 we show Spec2 {H; 60) for the first four and the 
last four pairs a,j in table [H respectively. We found the data for this 
and all the other pictures of second order spectra in this paper, by 
adapting in Matlab the Maple algorithm "reson25" included in the 
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appendix-a of the electronic version of [2 . Notice that there is always 
a point in the second order spectrum which approximates the isolated 
eigenvalue A^. We reproduce the coordinates of these points in the forth 
column of table H Nonetheless some of these values have imaginary 
part different from zero, in all cases the real part of the approximation 
of \a is accurate up to the 6*^ digit. This is confirmed by further 
numerical experiments. 

From figures |21 and 121 we can also say something about the shape of 
the second order spectra of H. Since V is of rank one, roughly speaking 
Spec2 {H; 60) should not be too far from the set 

{zeC : detTeoiP) = 0} 

which is shown in figure ^ (the matrix is tridiagonal and it is not 
very large so any of the standard computer packages can produce the 
data for this picture). According to corollary El the convergence of this 
sets to the interval [—2, 2] is o{k~'^) for all < g < 1 as /c — oo. 
Compare with figures El and 01 Notice that in the latter there are some 
perturbed points very close to the real axis inside the elliptical region 
which comprises most of the second order spectrum. Further numerical 
experiments confirm that they increase in number as j increases. Their 
real part are related to the non-real roots of q{x). These non-real roots 
of q{x) are generalized eigenvalues which according to Davies ^ can 
be regarded as (absolute) resonances of H. It would be interesting to 
explore further how these resonances are related to Spec2 {H\ k) . 

In figure jH we consider Spec2 {H; k) for a = 3 and j = 3 as /c in- 
creases. Notice that the convergence to the continuous spectrum seems 
to be much slower than the convergence to the eigenvalue. The point 
(pk £ Specg {H; K) which is furthest from Spec {H) seems to have real 
part approximately equal to 0. Figure El is a log- log graph of the imag- 
inary part of (pk when k varies from 100 to 1500. The slope of the line 
is very close to -1. This suggests that corollary El fails for g > 1. In 
table 121 the slope of the line in figure El changes from test point to test 
point. This suggests that the convergence to the spectrum might be of 
the order log"(fc)/fc for some a > 0. 
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k 


Slope between k and 




fc + 100. 


100 


-0.8494 


200 


-0.8590 


300 


-0.8647 


400 


-0.8688 


500 


-0.8719 


600 


-0.8744 


700 


-0.8765 


800 


-0.8782 


900 


-0.8798 


1000 


-0.8811 


1100 


-0.8823 


1200 


-0.8834 


1300 


-0.8844 


1400 


-0.8853 



Table 2. Slope between the test points k and k + 100 in figure El 
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Figure 1. Points where detT6o(/?) vanishes. 
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Figure 2. Spec2 {H; 60) selected values of a and j cor- 
responding to table n 
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Figure 3. Spec2 {H; 60) selected values of a and j cor- 
responding to table ^ 



0.08 



0.06 



0.04 - 



0.02 




-0.02 - 



-0.04 



-0.06 



-0.08 - 



-2 -1.5 -1 -0.5 0.5 1 1.5 2 2.5 3 

Figure 4. Spec2 {H; k) for a = 3, j = 3 and three 
values of k. 
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Parameter k 



Figure 5. Log- log plot of the maximum distance be- 
tween Specg {H; k) and Spec {H) for a = 3 and j = 3 as 
k increases. 



